Scattering and resonance of ultracold atoms with spin-orbit coupling 
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We derive an analytical expression for the scattering amplitude of two ultracold atoms of arbitrary 
spin and with general spin-orbit (SO) coupling, on the basis of our recent work (Phys. Rev. A 86, 
053608 (2012)). As an application, we demonstrate that SO coupling can induce scattering resonance 
in the case with finite scattering length. The same approach can be applied to calculate the two- 
body bound state of SO-coupled ultracold atoms. For the ultracold spin-1/2 fermi gass in three- or 
two- dimensional systems with SO coupling, we also carefully examine the renormalization relation 
of effective contact interaction with momentum cutoff, as well as the applicability of Huang- Yang 
pseudo-potential. 
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I. INTRODUCTION 

In the study of ultracold gases it is very important 
to understand the low-energy scattering properties of 
atoms. First, the two-body and three-body collisions are 
the underlying physics of many important experimental 
phenomena, e.g., the two-body decay and three-body re- 
combination. Second, understanding the behavior of the 
low-energy inter-atomic scattering amplitudes is indis- 
pensable in designing the effective inter-atomic interac- 
tions (e.g., the Huang- Yang pseudo-potential Bethe- 
Peierls boundary condition J2[ and the contact interac- 
tion with momentum cutoff [3[), which are widely used in 
the theoretical calculations. In the low-energy cases the 
effective interaction and the realistic interaction potential 
should lead to the same two-body scattering amplitude. 

In recent years, a class of synthetic gauge fields and 
spin-orbit (SO) coupling has been realized in ultracold 
Bose gases 0-QiJ and degenerate Fermi gases [ll], 03 
with Raman beams [ill 13 • A considerable amount of 
theoretical interest has been stimulated to understand 
the SO-coupling effect in both many-body [j"5H37j and 
few-body physics (37rl47j . 

Recently, the scattering amplitude of two ultracold 
fcrmonic atoms with Rashba SO coupling in the two- 
dimensional (2D) or quasi-2D confinements is calculated 
by two of us and coauthors [4l| . In addition, in Ref. [I?} 
we study the short-range behavior of the scattering wave 
functions of two SO-coupled ultracold atoms in a 3D uni- 
form system, and derive a modified Bcthc-Peierls bound- 
ary condition. Based on these results, in this paper we 
provide a systematic investigation on the two-body scat- 
tering amplitude of SO-couplcd ultracold atoms in 3D 
uniform space. We derive an analytical expression for 
the scattering amplitude of two atoms with arbitrary spin 
and SO coupling. Our approach can also be used to cal- 
culate the low-energy bound state of two spin-1/2 atoms 



with SO coupling. Furthermore, we show that the SO 
coupling can induce the scattering resonance. Namely, 
for the atoms with finite scattering length, the thresh- 
old scattering amplitude diverges when the SO-coupling 
intensity assumes some particular value. For the SO- 
couplcd spin-1/2 fcrmonic atoms, we also investigate the 
renormalization relation of 3D and pure-2D effective con- 
tact interaction with momentum cutoff, as well as the 
applicability of Huang- Yang pseudo-potential. We find 
that the form of the 3D and 2D renormalization relation 
is not changed by the SO coupling. Nevertheless, the 
Huang- Yang pseudo-potential cannot be directly used in 
the presence of SO-coupling. 

The remainder of this manuscript is organized as fol- 
lows. In Sec. II, we derive the exact analytical expres- 
sion for the scattering amplitude of SO-coupled ultracold 
atoms, and show the approach for the calculation of two- 
atom bound states. In Sec. Ill, we illustrate the scat- 
tering resonance induced by SO coupling. The renormal- 
ization relation of effective contact interaction and the 
applicability of Huang- Yang pseudo-potential is investi- 
gated in Sec. IV. The main results are summarized and 
discussed in Sec. V, while some details of our calculations 
are explained in the appendixes. 



II. SCATTERING AMPLITUDE AND BOUND 
STATE OF SO-COUPLED ATOMS 

A. spin-1/2 fermionic atoms 

We first consider the scattering amplitude of two SO- 
coupled spin-1/2 fermonic atoms. Without loss of gener- 
ality, the single-atom Hamiltonian of the system can be 
written as 



H 



P 2 



ii.) 



AM • P 



(1) 
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where P is the atomic momentum, and M and Z arc 
operators in spin space (here we have used h = 1 and the 



2 



atomic mass m = 1). The term AM • P describes the SO 
coupling and Z accounts for the residual spin-dependent 
part. Here, the eigen-values of M are of the order of 
unity, and A indicates the intensity of the SO coupling. 
For effective spin-1/2 systems in Refs. [U EH, one has 
A = 2k ri M = (<r z ,0,0) and Z = 5a z /2 + Q.a x /2, where 5 
is the two-photon detuning, k r is the recoil momentum, 
f2 is the Raman-coupling strength, and a z = \ t) (t | — | 4- 
) (I | and = | t) (4- | + | 4-)(t I are t ne Pauli operators, 
with |t) and |4) the spin eigen-states. 

The total Hamiltonian of the two atoms is given by 

H ib + H ib + U (f)i whcrG H ib (* = !' 2 ) is for thc 
ith atom, and U (r) is the spin-dependent interaction 

potential between the two atoms with relative position 

r = (x,y,z). Because the total momentum of the two 

atoms is conserved, the relative motion of these two can 

be separated from their mass-center motion. In the r- 

representation, the Hamiltonian for the relative motion 

is then 

H=f + Xc-p + B(K) + U(r) = H Q + U(r), (2) 

where p = — iV is the relative-momentum operator of the 
two atoms. The total momentum K = P^> + P^ of the 
two atoms is conserved during the scattering process, and 
behaves as a constant in our calculation. The operators 
cand B(K) read 



M (1) -M (2) , 



(3) 



B{K) = zW + zW + -K-(mW+mW). (4) 



In the stationary scattering theory, the incident state is 
regarded as the eigen-state of the Hamiltonian Ho for the 
free motion of the two fermionic atoms. With the Pauli's 
principle being taken into account, such an incident state 
for two spin— 1/2 fermonic atoms can be expressed with 
the spinor wave function 



|*f(r)) 



lk-r 



-ik-\ 



Air 3 /' 



\a,k) 



4^3/: 



-P 12 |a,-fc), (5) 



with k the relative momentum of the two atoms and P12 
the permutation operator of the spin of thc two atoms. 
Here the state \a, k) (a — 1, 2, 3, 4) of the two-atom spin 
is defined as the ath eigen-state of the operator ho(k) = 
Xc-k + B(K), with eigen-energy £(a, K ,k). The symbol 
satisfies a > a' when £(a,K,k) > £(a',K,k). In this 
paper, we denote 



t = (a,K,k) 



(6) 



as the set of these three quantum numbers. It is easy to 
prove that |^ t (r)) in Eq. ([5]) is an eigen-state of Hq with 
cigen-energy 



E t = k 2 +£(a,K,k). 



(7) 



In this paper, we assume that U(r) is a short-range po- 
tential with effective range . In the region r = \r\ > , 
we have U (f) ~ and further the low-energy scattering 

state \^[ + \f)} with respect to the incident state ^^(r)) 
can be expressed as [4l|, Hl| 



vJ/ t (O) (r-))+i? t Go(i?t;r,0)|S), (8) 



with |S) = (|t)i|4>2 - \i) i|t> 2) /\/2 the singlet spin state. 
Here B t is a r-independent constant and related with 
I'Jjj (r*)) through the relation 



dr'U(r')\^ + \f'))=B t \S). 



(9) 



For details, see Appendix A. In Eq. 
function Go (rj; r, f ') is defined as 



the free Green's 



G (v,r,f') = 



1 



rj + iO+-H 



5(f-f') 



(10) 



and is a (f, r ')-dependent operator for the two-atom 
spin. In this paper we consider the low-energy case with 
k <C 1/r*, and further assume that the SO coupling is 
weak enough so that A << Furthermore, we have 

proved [4?} that, in the short-range region 7% < r <C 1/k 

the function \^[ + (r)) behaves as 



|^ +) (r)) oc 



1 



\S)-i^S- 



|S). (11) 



Here, the scattering length a R is determined by both the 
detail of the potential U(r) and the SO coupling. In some 
special systems, e.g., the systems with \U(r), c| = and 
those in the current experiments [4HlO|, the scattering 
length a R is independent of the SO coupling and takes 
the same value as the scattering length in the systems 
without SO coupling |47j . 

We can now calculate the the coefficient B t in Eq. © 
and the inter-atomic scattering amplitude. As in 
Refs. HI HI], B t can be obtained from Eq. (JTTJ and 
the behavior of Go (rj;f,r') |S) in the short-range region 
< r << l/k. With calculations shown in Appendix 
B, we find that 

Gofer-, 0) « -l-^+i^+F^+i^c.^ 

(for r* <r « l/k), (12) 
where the operator F (77) is defined as 



H dk 



: „ ( \a",k")(a",k"\ \a"X){a",k"\ \ 
77 + i0+ - E t n rj + i0+ - k" 2 J ' 

(13) 
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with t" = (a",K : k"). Substituting Eq. ([12]) into Eq. flSJ}, 
we get the expression for ^^'(r)) in the short-range 
region 

A_ 



\¥ t +) (^) = \^ 0) (0))+iB t ^c- [ -) |S) 



-Bt 



1 /I 
-— \-+iE t 
4ir V r 



F(E t 



IS) 



(for r» < r << 1/fc) .(14) 

Comparing Eq. (fl"4"|) to Eq. (fTTj) and using the fact that 
F(?7)|S) <x |S), we obtain 



B, 



l/a B +iEl /2 -4ir(S\F(E t )\S) 



(15) 



According to scattering theory [49| , the scattering am- 
plitude f(t' <— t) between the incident state 
and an energy-conserved output state \^/[?\r)) with t' = 
(a' , K, k') is defined as 

f(t' f- t) = -2^ 2 / dr(*^(r)|C/(r)|^ +) (r)). (16) 



It is pointed out that, since U(r) ~ in the region 
r > r*, the integration in the right-hand side of the 
above equation is only done in the region r < . Under 
the low-energy condition k << 1/r*, in this region we 



have (*), J (r)| w (%,'(0)\, and then /(f <- t) can be 
re-expressed as 



(0), 



f(t'^t) = -2tt 2 <^ 
,(o) 



dr[/(r)|^ 0) (f)> 



27r^(0)|S)B t! 



(17) 



where we have used Eq. ([9]). Thus, using Eq. (|15j) . we 
finally have 



- (2tt) 3 <vI/W(0)|- 



l/o R + zS t 1/2 - 4 7 r J F(E t ) 



(0)) , 
(18) 



1 /2 

with i? ( = «y |-Et| for E 4 < 0. This is the exact analyt- 
ical expression for the low-energy scattering amplitude of 
two spin-1/2 fermonic atoms with SO coupling. 



where ± are for the systems of bosonic an fermonic 
atoms, respectively. 

The scattering wave function with respect to incident 
state \^[ 0) (r)) can still be denoted as \^[ +] (r)}. In the 
region with r > we have (Appendix A) 

|*| +) (r)) « |*fV))+Go(£t;r,O)|0) (20) 
with the r-independent spin state \<p) satisfies 



= / dr'U{r')\^{ + \r)) 



(21) 



This is very similar to Eq. ©, but the spin state \cj>) is not 
unique. Instead, \<p) can be different for different incident 



state (r)). Furthermore, as shown in Ref. [47[, in the 



short-range region r* < r << 1/fc the scattering wave 



function \^[ + \r)} behaves as 



* { t +) (?))=[^-A R )\x)-i'i 



(22) 



Here \x) is another r- independent spin state, and An is 
a r-independent operator in the spin space, which is also 
determined by the detail of the interaction potential U(r) 
and the SO coupling. For the cases of spin— 1/2 fermonic 
atoms, we have An — 1/an, and Eq. (f22j) reduces to 
Eq. (|llj). As in the above section, in the current experi- 
ments [■i-flOll for bosonic atoms with ID SO coupling, An 
is independent of the SO coupling. For instance, for the 
ultracold gases with spin-1 87 Rb atoms, we have 



A R = —V F=0 
a 



"2 



-V F = 



(23) 



where ao (0,2) is the scattering length with respect to the 
total atomic spin F = (F = 2) and Vf=o,2 are the 
associated projection operators. 

Using Eq. (fT2"j). we can obtain the expression of 

(r)) in the short-range region. Comparing such an 
expression with Eq. (f22j) , we have 



= -4tt|x) = 



4/T 



A K + iE} /2 -4nF(E t ) 



^ 0) (0)) (24) 



and 



B. Atoms with arbitrary spin 

The above approach can be straightforwardly general- 
ized to the general case of two fermonic or bosonic atoms 
with any kind of SO coupling and arbitrary spin. In these 
cases, the single-atom motion and the relative motion of 
the two atoms are still given by Eqs. (JIJ and ([2]), respec- 
tively. The incident state can be expressed as 



n ik-'. 



ik-r 



47r 3/2 



\a,k) ± 



47r 3/2 



Pi 2 |a, -k) 



(19) 



f{t' <- t) 



At. 



J . „,„J ^ 0) (Q)) 



iE t — AnF(Et) 



(25) 



Now we consider a simple case that An = 1/a with 

a a constant c-number and \a\ is much smaller than 

1/2 

the eigen- values of the operator iE t ' — 4irF(E t ). In 
this case Eq. (|2~51 can be simplified as f(t' 4— t) « 

-a(27r) 3 (*|! )) (0)|*f ) (0)). This approximate result is 



4 



also obtained in Ref. @ with Fermi's golden rule. On the 

other hand, if the eigen-values of is comparable to or 

1 li 

smaller than the ones of the operator iE t — AirF(Et), 
the contribution from that operator becomes significant 
and the Fermi's golden rule is no longer applicable. 

It is pointed out that, in the presence of SO-coupling, 
the state l^ "* (0)) depends on the quantum number 
t = (a,K,k), and thus changes with the direction of 
the atomic relative momentum k. Therefore, the scatter- 
ing amplitude f(t' «— t) in Eq. ([25]) is anisotropic with 
respect to the directions of the incident momentum k 
and output momentum k! . This anisotropicity is also 
observed in the experiment by I. B. Spielman et. al. @. 



200 




X (in units of k 1 ' 2 ) 



200 



C. two-atom bound state 

In the above subsections we derived the analytical ex- 
pression of the scattering amplitude of two ultra-cold 
atoms with SO-coupling. The Green's function approach 
used in our calculations can also be applied to derive 
the low-energy bound state of two SO-coupled atoms. 
We denote the energy of the bound state as Eb, then the 
binding energy can be defined as -Ebindmg = — (Eb — Eth), 
where E t h is the energy of scattering threshold or the low- 
est eigen-energy of Hq. As shown in Appendix A, when 
-^binding << r * 2 , the wave function \ty b (r)) of low-energy 
bound state is given by 

|* 6 (f)) « A 6 G o (£ 6 ;r,O)|0 b ), (26) 

in the region r > . Here Nb is the normalization factor 
and \<f>b) is a r- independent spin state. Furthermore, 



in the short-range region < 



« R 



-1/2 
binding 



1*6(0) 



has the same behavior as the low-energy scattering state 
l 1 !^ (r)) and thus can be expressed as |\E'b(r)) = — 
A R)\Xb) - i(X/2)c- (f/r)\xb), with \xb) a r-indepedent 
state in the spin space. As in the above section, with this 
fact and the short-range behavior of Go (rj; r, 0) given by 
Eq. (fT2|) . we find that \xb), \<Pb) and Eb can be obtained 
by 



-200 
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FIG. 1: (Color online) The threshold scattering scattering 
amplitude f(t 4— t) of two spin- 1/2 fermoinc atoms with ID 
SO coupling. We plot the f(t <- t) as a function of SO- 
coupling intensity A with the scattering length a R = 2k~ 1 '" 2 
(a) and a function a R with A = 1.25ft 1 ' 2 (b). Here we have 
used the natural unit with h = m = 1 and set the total 
momentum of the two atoms to be K = (Ko, 0, 0), with Ko = 
(black solid line), 0.3k 1 / 2 (red dashed line), 0.5k 1,/2 (green 
dotted line) and 0.8k 1//2 (blue dashed-dotted line). 



III. SO-COUPLING-INDUCED RESONANCE 



-iE\ 12 + AttF (E b ) 



\ X b) = A R \ Xb ), (27) 
|X6) = ~N b \4>„), (28) 



1 /2 I 

with E b ' = iy/\Eb\ior Eb < 0. In particular, for two 
spin— 1/2 fermonic atoms, we have = l/a R and \4>b) = 
|S). Then Eq. (J2TJ) becomes 



iEl /2 +ATT(S\F(E b )\S) 



(29) 



In above discussions, we derive the general analytical 
expression (|2"5| for the scattering amplitude of ultracold 
atoms with SO coupling. As an application, in this sec- 
tion we will show that SO coupling can induce scattering 
resonance of two ultracold atoms. For simplicity, here 
we only consider the scattering of two spin- 1/2 fermoinc 
atoms. 

In such a system, the incident state |^ (r)} takes the 

form in Eq. ©, and one has |^^ 0) (0)> oc |S). Thus, the 
inter-atomic scattering amplitude f(t' t) in Eq. (TTJ 
can be re-writtcn as 



With this equation one can obtain the bound-state en- 
ergy E b . 



5 




1/a R (in units of k 



FIG. 2: (Color online) The absolute value of scattering am- 
plitude of two spin-1/2 fermonic atoms with ID SO-coupling. 
Here we plot \f(t' < — t)\ for the cases with K = 0, A = 
1.25k 1/2 and E t > E th - We take t = (2, 0, 0) and t' = to (blue 
solid line), t\ (green dashed line), ti (red dotted line) and % 
(black dashed-dotted line), with to, 1,2, 3 defined in Sec. III. 



with the function d(A,E t ) defined as 



d(\,E t )=iEl /2 -47r<S|F(£7 t )|S>. 



(31) 



We first consider the case of threshold scattering with 



E 



In this case, we have 



t = E th . In this case, we have (f)) = |*| 0) 
and <i(A, E t h) usually takes a real value. Therefore, when 
the condition 



1 



d(A, E\ 



ih J 







(32) 



is satisfied, the threshold scattering amplitude diverges 
and a scattering resonance occurs. Since the scattering 
length a R can be finite, such a resonance is induced by 
the SO coupling. Comparing Eq. (|3~2"j) to Eq. (|2^|) . we 
further find that a bound state with zero binding energy 
appears at the resonance point. In experiments, one can 
observe the SO-coupling-induced resonance by tuning the 
scattering length or the SO-coupling intensity [50I ] . 

In Fig.[T]we plot the amplitude f(t <— t) of the thresh- 
old scattering of two spin-1/2 fermoinc atoms with ID 
SO coupling, i.e., M = (a z , 0, 0) and Z = Ka x , as in cur- 
rent experiments. In our calculation we assume the total 
momentum K of the two atoms is along the x-direction, 
i.e., K = (Kq,0, 0). The scattering amplitudes with re- 
spect to different values of Kq arc illustrated versus l/a R 
and the SO coupling intensity A. The appearance of res- 
onance is clearly shown. 

Next, when the scattering energy E t is larger than 
the scattering threshold, the function d(X,E t ) takes a 
complex value. Thus, the scattering amplitude fit! <— 
t) cannot divergent. Nevertheless, as a function of 
l/a R , the absolute value of the scattering amplitude 



still achieves a local maximum value when the condi- 
tion Re[l/a R + d{X,E t )} = or l/a R = -Re[d(A, E t )\ 
is satisfied. To illustrate this effect, we also calculate 
the scattering amplitude f(t' 4— t) for two spin-1/2 fer- 
moinc atoms with M = (cr z ,0, 0) and Z = na x . In 
Fig. [5] we plot \f(t' <— t)| versus l/a R for the cases 
with K = 0,t= (2,0,0) and t' = in,i,2,3 where t = t, 
h = (4,0, sfW ^?/\), t 2 = (4,0,VA 4 -k 2 /A-0.5k 1 / 2 ) 
andi 3 = (4,0, VA 4 - k 2 /A-0.8k 1/2 ). The peak behavior 
of the scattering amplitude is clearly illustrated. 

We mention again that the resonance discussed here 
is essentially induced by the SO-coupling term in the 
Hamiltonian Ho defined in Eq. ([2]) . From Eq. ([2]) one can 
easily find that such a term can be omitted when k >> A, 
and therefore the SO-coupling-induced resonance is sig- 
nificant only when the atomic relative momentum k is 
small enough i.e., k < A. Finally, we emphasize that the 
calculations in this seciton can be directly generalized to 
the systems of atoms with arbitrary spin and SO cou- 
pling, and the SO-coupling-induced resonance can also 
appear. 



IV. EFFECTIVE INTERACTIONS 

In the many-body theory of ultracold gases, the inter- 
atomic interaction is usually modeled by some simple 
effective potentials. The most widely-used effective in- 
teractions include the Huang- Yang pseudo potential [l[, 
the Bethe-Peierls boundary condition J3 and the contact 
interaction with a momentum cutoff [3j . 

For a given system, the inter-atomic scattering ampli- 
tude given by the effective inter-atomic interaction should 
be the same as the one from the realistic interaction po- 
tential U(f). In our previous works [U [47|, we have 
shown that to satisfy such a condition, the Bethe-Peierls 
boundary condition for 3D ultracold g well as 

the renormalization relation for the contact interaction 
of quasi-2D gases, should be modified in the presence of 
SO coupling. In this section, we consider the contact in- 
teraction with a momentum cutoff in a 3D and pure-2D 
uniform system of the SO-couplcd spin-1/2 Fermi g cLS, cLS 
well as the Huang- Yang pseudo potential. 



A. contact interaction in 3D system 

In a 3D system, the contact interaction with a momen- 
tum cutoff takes the form 



,t 



Uo 

£3 2-4 a q/2+k,t U 'q/2~k,l" , q/^-k',i" 1 q/2+k',f 
k.k' .q 



(33) 

Here, at^ and ap i(7 are the creation and annihilation op- 
erators for an atom with momentum p and spin a, L is 
the length of the system, and the summation y^g g, - is 



6 



done under the condition max(|fc|, \ k'\) < k c with k c is a 
cut-off momentum. 

The rcnormalization relation for this contact potential, 
i.e., the relationship between Uq and k c , can be obtained 
from the condition that U e g and the realistic inter-atomic 
interaction should lead to the same low-energy scatter- 
ing amplitude. For the cases without SO-coupling, the 
standard calculation gives the well-known result 



1 



1 



1 



(2tt) s 



dk"- 



1 



(34) 



|fc"|<fe c 



with a s the s-wave scattering length. 

In the presence of SO coupling, the correct renormal- 
ization relation can be obtained by the same procedure. 
The effective scattering amplitude f c s(t' «— t) given 
by U ft can be obtained from the Lippmman-Schwingcr 
equation for the two-body T-operator T e g (77) with re- 
spect to U e g: 



T c s (17) = U eS + U cS G (77) T ff (v) > 
where the Green's operator Go (77) is defined as 

1 



Go (77) 



T] + iO+-H 



(35) 



(36) 



Now wc derive the expression of the scattering ampli- 
tude f e ff(t' <— t) given by the effective interaction U e g. 



To this end, we introduce the notation 

1 



1*0 



ik-r 



(27T) 



3/2 



(37) 



for the plane wave in the spatial space. We also define 
inner product for the wave functions in the spatial space 
as 



(k'\k) = S(k-k' 



(38) 



With these notations, the first-quantization form of [7 e ff 
and Go (77) can be re-expressed as 



U, 



off 



Up 

(2tt)'° J\k\,\k>\ <kl 



\k)(k'\®\S)(S\dkdk', (39) 



and 



Go (v) = £ J * 



|fc)(fc| ® \a,k)(a,k\ 
77 + «0+ - E t ' 



respectively. Here we have t = (a, K, k) as before. 

Using the Lippmann-Schwinger equation, we can ob- 
tain the equation for the T- matrix clement: 



T eff ( v ,t',t) = 



Up 
(2tt) J 



{a',k'\S)(S\a,k) 



Up 
(2tt) ; 



dk"~ 

|fc"|<fc c V 



(S\a",k") 



E t 



-T cff (r),t",t) 



(40) 



where we have t' = (a',K',k') and T e g (77, t' , t) = 
(k'\(a',k'\f eS (77) \a,k)\k). The notations t" and 
T e g (77, t", t) are defined similarly. From Eq. (|40|). we find 
that T e g (77, t' , t) can be expressed as 

T eS (r,,t',t) = J^-( a ',k'\S)u (f],a,k) (41) 

(Z7r) ^ ' 



where 7/(77, a : is independent of 0/ and k' . Substituting 
Eq. gl]) into Eq. (|40]), we get the result 



f 77, a, kj 



(S\a,k) 



1 ~ jBf (So" Jj 



dfc' 



(S|a",fc"> 



»" J \k" <fe c " 77+10+-B 



Substituting Eq. (|4"2")l into Eq. (|4*Tj) . we can obtain the 
expression of the T-matrix element T e g (77, t' , t). 



The scattering amplitude f c e(t' «— i) for two spin- 1/2 
fermionic atoms is defined as 



fes(t' <- 1) = -27T 2 (W\T eS (E t ) |*H>, (43) 

where the incident state is }) = 2~ 1/2 [|a, fc)|fc) - 
(Pi2|a, — k))\ — k)], with P12 the permutation operator 
for the spin of the two atoms. Therefore, it is apparent 
that we have 

fes(t' <- t) = -2tt 2 [T cff (E t ,t', t) - T ' ff (E t , t', t)} (44) 

with T' eS (E u t',t) = (k'\(c/,k'\f eS (r,)[P 12 \a,k)}\ - k). 
Substituting our result of T e g (E t ,t',t) into Eq. (|%4"|) . wc 
have 



(42) f eS (t' <- t) = 



%(a>,k'\S)(S\a,k) 



1 _ JJo I v f,- , ,r„ |(s|q".fc")r 

1 (5^ I ^a" J|fc''|<fc c " K ,7+i0+-£; t ,, 



(45) 



7 



By requiring that the zero-energy effective scattering 
amplitude f e g(t' <— t) be equal to the realistic scattering 
amplitude, i.e., 



/off <" t) = f(f <- t) 



(46) 



with E t = and t) given by Eq. ([15]). we obtain 

the following renormalization relation for systems with 
SO coupling in the limit k c — > oo: 



1 



1 



1 



47ra R U (2tt) 3 Jft,^ \k"\ 



dk"- 



1 



(47) 



coupling. The straightforward calculation yields 
(t> <- t) 

4^(*i 0) (0)|S)(S|^ 0) (0)> 



1/2 



C - In ^ - 27rJ"( 2£l ) (£ t ) 



,(50) 



0.5772... is the Euler T number, £ t is the 
^ 2D ^ with respect to the eigen-state 



where C 

eigen-value of Hq 
\^[°\f)). In Eq. ([HI the function (77) is defined as 

1 „ r - 2 



(2D) 



(V) 



(27T)" 



J dk"\(S\a",k" 



Here we have used Eq. (|B3j) . Comparing Eq. (|34|) with 
Eq. (|47|) . we find that for the 3D contanct potential with 
momentum cutoff, the form of the renormalization rela- 
tion is not changed by the SO coupling. In the presnece 
of SO coupling one only needs to replace the scattering 
length a s with a R . 

We mention that Eq. P?)) for the bound-state energy 
Ef, can also be obtained from the contact potential U e g 
with renormalization relation (|47p . 



B. contact potential in pure-2D system 

Our above discussion can be directly generalized to 
pure-2D ultracold gases of spin- 1/2 fermonic atoms with 
SO coupling. We assume the atoms are moving in the x— 
y plane. Thus, the single-atom Hamiltonian is also given 
by Eq. ([IJ, with P = (P x ,P y ) the single-atom momentum 
in the x — y plane. The Hamiltonian for the relative 
motion is then 

H (2D) = f 2 + m+ + U2o{ff) = H {2D) + UiD ^ 

(48) 

where K is the two-atom total momentum, p— (x,y) is 
the two-atom relative position in the x—y plane, and U2D 
is the two-atom interaction potential in the 2D space, 
with effecitve range p*. Here the operators c and B(K) 
in the two-atom spin space are defined as in Sec. II, 
and SO-coupling intensity A is also assumed to be small 
enough so that the condition A < < 1/p* is satisfied. 

Similar as in Sec. IIA, the incident state in the scat- 
tering, or the eigen-state of the Hamiltonian H^ D \ can 
be expressed as 



*i 0) (r)): 



237% 



\a, k) 



-ik-p 



2 3 / 2 7T 



PuK-A) 



(49) 



where the state \a, k) for the two-atom spin is defined 
as in Sec. IIA, and we have t = (a, K, k) as before. In 
Ref. [4l| , we have calculated the 2D scattering amplitude 



f 2D \t' <- t) between |*J u, (f))and |#J, ; (r))for the cases 
with Rashba SO-coupling. The method applied there can 
be directly used for the cases with arbitrary type of SO- 



V 



Z0+ - Ef 



n 



i0+ 



(51) 



with t" = (a", K, k"). The parameter cLr in Eq. ([50)) can 
be determined by the following condition (sec Appendix 
B of Ref. [4l| where dp is denoted by d): in the region 
with < \p\ << 1/fc, the solution \ipn(p)) 01 equation 

T{p)H {2D) T^{p)] \MP)) = E t \Mp)) (52) 

satisfies \iPt(p)) (In \p\ — ln<iT)|S). Here the rotation 
T{p) is defined as T(p) = exp{i\c x x/2) exp(iAc y y/2). 

Now we consider the 2D contact potential which takes 
the form 



U, 



(2D) 



u, 



(2D) 



off 



L 2 



^ a. 



,t 



g72+fc,t a 972-fe4 a 972-fc'^ a <f/2+fc',f 



Here, and 



(53) 

are the creation and annihilation op- 
erators for an atom with momentum p and spin a . L is 
the length of the system, and the summation £, ^ is 

done under the condition max(|/c|, |fc'|) < k c with k c is a 
cut-off momentum. 

As in the 3D case, the renormalization relation for 
this contact potential can be obtained from the condi- 
tion that f7 G ^- D ' and the realistic inter-atomic interaction 
should lead to the same low-energy scattering amplitude. 
In the absence of SO-coupling, the standard calculation 
gives [HH 



1 



U, 



(2D) 



(2tt)' 



dk"- 

|fc"|<fc c £■ 



1 



\k"\ 2 ' 



(54) 



where the physical parameter e > is the binding energy 
of the 2D two-atom bound state for the cases without SO 
coupling. 

In the presence of SO coupling, the effective scatter- 
ing amplitude f e s(~tf <— t) given by U c g can be obtained 
with the method used in Sec. IV. A. The straightforward 
caclulation gives 



A2D) , , 
Jeff l l 



t) 



4 7 r 3 (^ ) (0)|S)(S|vl>r(0)) 



(0), 



1/2 



ln£; c - 27rJ"(2D) (£ t ) 



(55) 
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By requiring that f^ D) {f <- t) = f& D \t' <- t) with 
f( 2D \t' <— t) given by Eq. (|50l) . we obtain the following 
renormalization relation for systems with SO coupling in 
the limit k r — > oo: 



U, 



(2D) 



dk" 



1 



|fe"|<fec 



SR+\k"\2' 



(56) 



with en = 4exp(— 2C)/d 2 R . Comparing Eq. (|5"4"]) with 



Eq. (|56|) . we find that for the 2D contanct potential with 
momentum cutoff, the form of the renormalization rela- 
tion is not changed by the SO coupling. In the presnece 
of SO coupling one only needs to replace the scattering 
length e with sr. 



C. Huang- Yang pseudo-potential 

In the end of this section consider the Huang- 
Yang pseudo-potential Q for the SO-coupled spin-1/2 
fermi gas in a 3D system. The Huang- Yang pseudo- 
potential [l| is given by 



d 

U H y = 4na5 (r) — (r-) 
or 



(57) 



with a the scattering length. Since Uhy is a zero-range 
potential, the scattering state with respect to Uhy, i-e., 
the solution of the equation 

[H + U HY ] if) = E t \4 +) if) (58) 
with out-going boundary condition, should take the form 

\4 + Hf) = \^ Q) (^)+B t G (E t ;r,0) |S) (59) 

in the region with r = \r\ > 0, and thus satisfy the 
Bethe-Peierls boundary condition. In the absence of SO- 
coupling, the Bethe-Peierls boundary condition is 



lim|^ +) (f))oc (---) 



(60) 



According to this condition, the function r\ip[ + \r)) is 
continuous in the point r = 0. Therefore, the partial 
derivative d/dr in Uhy is well-defined for the function 

r\i>\ ; (r)}- Namely, the operation of the Huang- Yang Appendix A: Behavior of scattering wave function 

pseudo-potential Uhy on the wave function (r)) is 
well-defined. 

Nevertheless, in the presence of the SO-coupling, as 
shown in Ref. [47| , the Bethe-Peierls boundary condition 
is modified to 

hm |^ (+) (r)) «(---) |S)-4c- (-) |S). (61) 

Due to the anisotropic term — iXc ■ r/(2r)\S), the func- 
tion r\Tjj[ + \r)) is not continuous in the point r — 0. As 



a result, the partial derivative d/dr in Uhy is not well- 
defined for the function r\ipl + \r)), and thus the opera- 
tion of Uhy on the wave function \ipl + \r)) is no longer 
well-defined. Therefore, in the presentee of SO-coupling, 
the Huang- Yang pseudo-potential is not consistent with 
the modified Bethe-Peierls boundary condition, and thus 
cannot be directly used in the theoretical calculations. 



V. DISCUSSION 

In this paper we derive the analytical expression of 
the scattering amplitude of two ultra-cold atoms with 
SO coupling. Our approach can also be used to calcu- 
late the wave function and the energy of two-body bound 
state. Moreover, we show that the SO coupling can in- 
duce inter-atomic scattering resonance. The influence of 
such a resonance in the many-body physics of the SO- 
coupled ultra-cold gases remains to be explored. With 
the expression of scattering amplitude, we further prove 
that the renormalization relations of the 3D and 2D con- 
tact potentials with momentum cutoff are not changed by 
the SO coupling. It is pointed out that, this result pro- 
vides a solid basis for some previous theoretical works 
(e.g., Refs. HHHHHil) where such a renormalization re- 
lation is used for the SO-coupled gases without proof. 
We also show that in the presence of the SO-coupling, 
the Huang- Yang pseudo-potential is no longer consistent 
with the modified Bethe-Peierls boundary condition, and 
thus cannot be directly used in the theoretical calcula- 
tions. 
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In this Appendix we give the explicit definition of 
the effective range r* for the potential U(r), and prove 
Eqs. (0 HI HOI [2D EH) for the behavior of the scatter- 
ing state f)) and bound state |4 , f,(r)) in the region 
r > TV Although this problem has been discussed by 
us in Refs. J4?| and 41 1, the treatment there is not fully 
rigorous. Here we provide a more explicit analysis. 

We first consider the case of spin— 1/2 atoms. Accord- 
ing to the scattering theory, |^j + ' l (r)) is determined by 
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the Lippmman-Schwingcr equation 



^ +) (r)> = l*r«} + 



(0), 



df'G (E t ,r,r')U(r')\^ + \r')), 

(Al) 

with the Green's function Go defined in Eq. (JTUJ). Here 
we assume the potential U(r) is negligible when the inter- 
atomic distance is larger than a characteristic length ro, 
i.e., U(f) ~ in the region r > r . Thus, the integration 
in Eq. (|A1[) is only effective in the region r' < r . 

If the momentum fc is low enough so that k << 1/Vo, 
when r — > oo and r' < ro, the function Gq(E, r, r ') varies 
very slowly with respect to r' and we have Go{E, r, r ') rj 
Go(E,r,0). Therefore, in the limit r — > oo, the solution 
of Eq. (|A1[) takes the form 

|^ +) (r)) « \* { t °\r))+G Q (E u r,0M, (A2) 

where the spin state |0) is related to |^j + ^(r)) via the 
equation 



= / dr'U(r')\y { t + \r')). 



(A3) 



the condition ro < < 1/fc is satisfied and the anistropicity 
of U(r) and Un(r) is small enough so that the relevant 
effective ranges of U and U-r. are in the same order of 
magnitude with ro- 

Our above analysis can be directly generalized to the 
general cases of atoms with arbitrary spin. In that case 
the scattering state \^[ + \r)) takes the form in Eq. (|A2j) 
in the whole region with r > r*, and the r-independent 

state \(f>) is also related to (r)) via Eq. (|A3j) . Nev- 
ertheless, now \4>) is not unique. It can be different for 

different incident states \^[°\r)). Then we can obtain 
the result in Eqs. pD l ET jl . 

Now we consider the low-energy bound state of two 
SO-coupled atoms. It is clear that the wave function 
\^b(r)) of the bound state is given by the equation 

\V b (r)) =N b [ dr'G (E b ,f,r')U(r')\* b (r')) (A7) 



with E b the energy of the bound state and N b the normal- 
ization factor. With similar discussion as above, when 
|-Ebinding| << we have 



Furthermore, due to the facts P 12^1 (—r)) = 

-\$>\ + \r)) and Pi 2 U(-r)P 12 = U(r) with P 12 the per- 
mutation operator of the spin of the two atoms, one finds 
that P 12 [U(-r)\^ { f +) (-r))} = -U(r)\$[ +) (r)). This re- 
sult yields 



\<f b (r)) ~ N b G o (E b ,r,0)\<i>b) 



(A8) 



in the region r > r*. That is the result in Eq. ([25)1 . In 
this paper we assume the condition |i?binding| << 1/?"* is 
satisfied. 



df'U(f')\^ + \r')} = Bt\S) , (A4) 



with B t a constant number. Therefore, when the inter- 
atomic distance r is large enough, we have 



|*<°>(f)) 



B t G o (E u r,0)\S). 



(A5) 



A direct result is that there exists a characteristic 
length r* and Eq. (|A5I) is applicable in the region r > r* . 
It is apparent that we have r* > ro. In this paper we de- 
fine such a characteristic length as the effective range of 
the potential U{r). Therefore, Eqs. (SI EH) or Eqs. (SJ 
[9]) are naturally satisfied when r > r* . 

When U (r) is spherical, in the low-energy case one only 
needs to take into account the s-wave scattering. The 
straightforward calculaiton shows that, as a result of the 
s-wave approximation, we have r* = ro. When U(f) is 
anisotropic, the effective range r* is larger than ro. In 
our current paper and Ref. [47] . we assume the kinetic 
energy of atomic relative motion is low enough so that 
the effective ranges of U (r) and the "rotated potential" 



U R (r) = 

a i\Ca:x/2 i\c y y/2 i\c z z/2i 



e ^a,-,/ „„, _g ^f/(y^ e — - »-/- e VI ~ e 



-iXc x x /2 — i\c y y /2 — i\c z z/2 



(A6) 



is much smaller than This condition can be satisfied 
in the dilute gases. In particular, in the systems where 



Appendix B: The short-range behaviors of free 
Green's function 



In this Appendix we prove Eq. (|12p for the short-range 
behaviors of Go(E t , r, 0) with the following three steps. 
First, with the fact 



it is easy to show that 



G (r),r,0)=J2 / dk 



,e^"- f \a",k' / )(a,k // \ 
(2tt) 3 r) + i0+ - E t » 



(Bl) 



(B2) 



with t" = (a" ,K,k"). Eq. (jB2j) and the completeness 
relationship J2 a " \ a " i k") (a" , k"\ = 1 lead to the result 



G o (r?,r,0) = / dk 



(2tt) 3 t] + i0+ - fc" 2 



%k" -f 



+ I dk"- ~\a",k")(a",k"\ x 

(2tt) 3 



T) + i0+ - E t „ r) + i0+ - fc" 2 



(B3) 



10 



It is pointed out that, in the limit r — > 0, the integration 
in the right-hand side of Eq. (|B3[) does not diverge. On 
the other hand, we also have 



ik" -r 



dk"- 



1 



(27r) 3 V + i0+ - k" 2 Anr 
with argr/ G (— it, ir]. Due to these facts, we have 

1 



lim G (?7, r, 0) 



+ 0{r u ). 



(B4) 



(B5) 



Second, as in Ref. (47]], we introduce a unitary trans- 
formation TZ(r) as 



7\L(r) = e iAc xK/2 e iAc„y/2 e iAc 2 z/2^ 



(B6) 



with c = {c x ,c y ,c z ). The rotated free Hamiltonian 
H R =1l{r) H TV{r) can be calculated as 

H m =f- 2Xd(Xr) ■ p + W(f) (B7) 
with operators d = (d x -, dy,d z ) and W given by 
d x (Xr) 
d y (Xf) 

d z (Xr) 



0; 

2 

c JL-n{r) C -^ll\r), 



K(r)fll\r); (B9) 



(B8) 
;(B9) 

(BIO) 



and 



W(r) 



V • d(Xr) + H(r)B(K)1l\r) + 
|d(Ar)| 2 -^(r)M!^(r)|. (Bll) 



Now we define the rotated Green's function 
G(m.(rj, f, 0) as the Green's function with respect to 
Hqr, i-e., 



G r(?7, r,0) 



1 



r? + i0+-F R 



<5(r) 



(B12) 



According to Eqs. (|B8tiB10|) . we have d(Xr) = 0(Xr). 
Namely, the SO-coupling term in Hr vanishes in the limit 
r — > 0. Due to this fact, in such a limit Gor(t], r, 0) has 
the same behavior with the Green's function Gw (?7, 0)) 
which is defined as 



Gw(?7,f, 0) = 



rj + i0+- [p 2 + W(0)\ 



6(r) 



oiyjri— W n r 

4irr 



(B13) 



with W„ and |W„) the n-th eigen- value and eigen- 
state of the operator VT(0), respectively. Since 
lim r _j.o Gw(rj, r, 0) is the sum of — l/(4-7rr) and a r- 
independent operator, we have 



lim G R (57,^,0) = \-G' Q (rj), 

r— >o 47rr 



with G' (t]) a r-independent operator in the space of two- 
atom spin. 

As in Ref. [13] , the result in Eq. (|B14|) can be proved as 
follows. According to the definition of Gor, it is obvious 
that 



G o (r,,f,0) =K(7)G OR (r),f,0). 



(B15) 



Due to this relation and Eq. (|B5j) . Gor(t?, r, 0) can be ex- 
pressed as the sum of —l/(4irr) and another term which 
converges in the limit r — > 0. Then we have 

oo 

G R(r ? ,f,0) = -— - + J2r n G' n (r,) 



n=0 



+E E E rn ^.™«(^) J| ™fa)' ( B16 ) 

1 = 1 mi=—l n=0 

with (r, 0, 0) the spherical coordinate. Here Y^ mi (9, 4>) 
are the spherical harmonic functions, G' n (rj) and 
Ji,mi.n(jf) are r-independent operators in the space of 
two-atom spin. 

In the region r > 0, the Green's function Gor(^, r, 0) 
satisfies 



- 2Ad(Ar) ■ p + W(r) G 0K (ri, r, 0) = ijGtafa, r, 0) . 



(B17) 



Substituting Eq. (|bT6|) into Eq. (jBT7| and comparing 
the coefficients of the term r~ 2 in both sides, we find 
that because d(Xr) = O(Xr), one has J;, TOi ,o(??) = 0. 
Therefore, in the limit r — > 0, Gor(?J, r, 0) behaves as 
in Eq. (|B"T4) . 

Third, using Eq. (|B15[) and Eq. (jB14j) . we finally have 

limGo^r-O) = - J- +i A«?. (-) +G^). (B18) 
r->o 47rr S7r \r J 

Now we derive the r-independent operator G'^rf) in the 
space of two- atom spin. According to Eq. (|B18j) , we find 
that 



G' (v) 



lim I G (ry,r,0) + -^ 



+ lin^G„(,,-r,0) + i i-) 



Using Eqs. (jB3() and (|B4) and the fact that 



lim 

r->0 



dkf(k)e lk ^ + J dkf{k)e- lk F 
we find that 

1/2 



(B19) 



dkf(k), 
(B20) 

(B21) 



(B14) 

with the operator F (rf) defined in Eq. ((T3) . 
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In the short-range region, the behavior of Go (77; f, 0) is That is the result in Eq. (fT"2|) . 
given by lim r ^o Go(rj, f, 0). Thus, substituting Eq. (|B2ip 
into Eq. ([B18]) . we obtain 

Gofaf.O) « -^(j:+V /2 )+F(,) +l Ae. (?) 

(for r* <r« 1/fc). (B22) 
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